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We determine the energy-level shift experienced by a neutral atom due the quantum electromag- 
netic interaction with a layered dielectric body. We use the technique of normal-mode expansion to 
quantize the electromagnetic field in the presence of a layered, non-dispersive and non-absorptive 
dielectric. We explicitly calculate the equal-time commutation relations between the electric field 
and vector potential operators. We show that the commutator can be expressed in terms of a gen- 
eralized transverse delta-function and that this is a consequence of using the generalized Coulomb 
gauge to quantize the electromagnetic field. These mathematical tools turn out to be very helpful in 
the calculation of the energy-level shift of the atom, which can be in its ground state or excited. The 
results for the shift are then analysed asymptotically in various regions of the system's parameter 
space - with a view to providing quick estimates of the influence of a single dielectric layer on the 
Casimir-Polder interaction between an atom and a dielectric half-space. We also investigate the 
impact of resonances between the wavelength of the atomic transition and the thickness of the top 
layer. 
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I. INTRODUCTION 

The question of the interaction between a neutral atom 
and a macroscopic dielectric body, once of purely aca- 
demic interest, has recently been promoted to a real- 
life physics problem thanks to the rapid developments 
in nanotechnology and experimental techniques. It is 
no longer the case that this interaction, the so-called 
Casimir-Polder interaction, is a tiny effect that can be 
ignored in all practical situations. Instead, on the length- 
scales that nanotechnology nowadays operates in, disper- 
sion forces, as they are also called, become significant and 
may appreciably influence miniaturized physical systems. 
Many of the current ambitions of cold-atom physics to- 
wards quantum computation and a variety of nanotech- 
nological applications involves the trapping and accurate 
guiding of single atoms above dielectric substrates, so- 
called atom chips. With these the nearby environment of 
a trapped atom usually consists of a complicated array 
of inhomogeneous dielectrics. The questions then aris- 
ing are: what are the magnitudes of the Casimir-Polder 
forces felt by the atom, and can one possibly engineer 
the types and shapes of surrounding materials either to 
minimize unwanted dispersion forces or to make them op- 
timally contribute to the trapping or guiding? In order to 
investigate such possibilities one needs to go beyond sim- 
ple featureless geometries and ground-state atoms and 
gain flexibility. The perhaps least sophisticated but still 
interesting example to study in this context is to consider 
a neutral atom, possibly excited, above a layered dielec- 
tric half-space, cf. Fig. Q] If the atom is in its ground 
state, then the Casimir-Polder force is always attractive 
for material surfaces with refractive indices greater than 
1. In such case it is desirable to derive simple analyti- 
cal formulae that would allow one to obtain quick esti- 
mates of the magnitudes of the forces involved in terms 
of the optical properties of the layer and the substrate 



On the other hand, if the atom is in its excited state, 
then, as it is widely recognised 0] , the potential acquires 
a oscillatory contribution that can result in a repulsive 
force. Additionally, the presence of the layer creates the 
possibility of a resonance between the wavelength of the 
atomic transition and the thickness of the layer, which 
could lead to a suppression or enhancement of the inter- 
action. 

There exist a variety of theoretical approaches devised 
to study the Casimir-Polder interaction (see e.g. Q for 
a recent list of references) but perhaps the most success- 
ful ones being the linear response theory [4| and phe- 
nomenological macroscopic QED [5j. By using linear re- 
sponse theory [4( and expressing the field susceptibili- 
ties in terms of Fresnel reflection coefficients @, Q , one 
can express the Casimir-Polder interaction as an inte- 
gral along the imaginary frequency axis of the product of 
the atomic and field susceptibilities. Thus in practice the 
problem is reduced to the calculation of the classical elec- 
tromagnetic Green's tensor expressed in terms of Fresnel 
coefficients. Such calculations, while straightforward in 
principle, tend to be quite tedious and often inevitably 
lead to the use of numerical methods. However, there 
is a benefit to studying problems in quantum electrody- 
namics by using physically transparent methods that do 
not obscure the basic underlying physics. For the kind 
of geometry of plane layered dielectrics considered in this 
paper, the technique of electromagnetic field quantization 
based on a normal-mode expansion 0] seems to be best 
emphasizing the physics of the problem, namely the fact 
that the system supports two kinds of modes of the elec- 
tromagnetic field [8|: these are travelling modes with a 
continuous spectrum and trapped modes with a discrete 
spectrum, i.e. occurring at only certain allowed frequen- 
cies. The trapped modes arise because of repeated total 
internal reflections within the top layer of higher refrac- 
tive index than the substrate, and emerge as evanescent 
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waves outside the wave-guide. This gives rise to an in- 
tricate assortment of evanescent modes outside a layered 
dielectric where evanescent waves with continuous spec- 
trum, also arising in a half-space geometry Q, are super- 
posed with discrete evanescent modes that arise only in 
the presence of the slab- like waveguide pj . In the frame- 
work we apply in this work, in the same spirit as e.g. 
0, Q, the use of standard perturbation theory renders 
all calculations explicit and it is possible from the outset 
to track down and remove if necessary any ambiguities 
that tend to remain hidden in more elaborate theories. 
For example, linear-response theory results in an integral 
over the Fresnel reflection coefficients but gives no indi- 
cation of whether the evanescent waves associated with 
the trapped modes contribute to the Casimir-Polder in- 
teraction or not. The question is answered at once if the 
normal-modes approach is used instead, see Also, 
interpretations of more complicated field-theoretical ap- 
proaches 10] can be put to an explicit test pj. 

The purpose of this paper is twofold. Firstly, it aims 
to support current experimental efforts by providing a 
range of analytical formulae useful for quick estimates 
of the dispersion forces acting on an atom placed in the 
vicinity of the layered dielectric, with particular emphasis 
on the corrections caused by the layer as compared to the 
standard half-space results reported in @. It also inves- 
tigates the resonant interaction between an excited atom 
and a layer in the search for the possible enhancement 
or suppression of the Casimir-Polder force. Secondly, it 
formulates a simple and explicit theory based on well 
understood concepts of theoretical physics such as per- 
turbation theory and electromagnetic field quantization 
in terms of a normal-modes expansion. The theoretical 
aspect, although serving only as a means to a practical 
end result, turns out to be interesting in its own right. 
The perturbative approach used in this work leads to the 
problem of the summation over the modes of the electro- 
magnetic field, which is non-trivial because of the dual 
character of the modes of the electromagnetic field. The 
task of adding the discrete and continuous field modes is 
elegantly accomplished by the use of complex-integration 
techniques. This allows us to explicitly show that the 
canonical commutation relations between the field op- 
erators are satisfied, which is equivalent to saying that 
the completeness relation of the normal-modes holds in 
the geometry considered. Although this is not a surprise 
because the field modes are solutions of a Hermitean op- 
erator's eigenvalue problem, the explicit calculation we 
carry out provides us with the mathematics necessary to 
complete a typical perturbative calculation in this geom- 
etry. It also allows us to cast the end result in a simple 
and elegant form that is easy to study analytically in 
various asymptotic regimes. The same technique could 
be applied to any similar perturbative problem is such a 
geometry. 

This paper is organised as follows. First we quantize 
the electromagnetic field in the presence of a layered di- 
electric, Section [TTJ Then, in Section III C[ we explicitly 



prove the completeness relation for the electromagnetic 
field modes. Equipped with the necessary mathematical 
tools, we proceed to calculate the energy shift in Sec- 
tion [Till and then study it analytically (Section IIV[) and 
numerically (Section [Vj) . 



II. FIELD QUANTIZATION IN THE 
PRESENCE OF A LAYERED BOUNDARY 
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FIG. 1: Atomic dipole moment in the vicinity of the lay- 
ered dielectric. The dielectric function is a piecewise constant 
function of the coordinate z. 



Our ultimate aim is to work out the energy-level shift 
in an atom caused by the presence of a layered dielectric. 
In order to obtain a result that fully takes into account 
retardation effects, the quantization of the electromag- 
netic field is necessary. To emphasize the physics of the 
problem we choose to quantize the electromagnetic field 
by a normal- mode expansion as described in [12]. The 
dielectric environment we consider (cf. Fig. [IJ consists 
of a substrate, a dielectric half-space occupying the re- 
gion of space z < —L/2 described by a dielectric constant 
e s = rig , and on top of that substrate an additional dielec- 
tric layer of thickness L, which has a dielectric constant 
e\ = nf. We assume that the dielectric constant of the 
layer is higher than that of the substrate e\ > e s in order 
to account for modes that are trapped inside the layer. 
Although we work with this assumption, the final result 
will turn out to be valid even when the reflectivity of 
the substrate exceeds that of the layer, but that is the 
physically less interesting case. Throughout this paper 
we shall assume all dielectric constants to be frequency 
independent so that the optical properties of the system 
are described solely by a pair of real numbers, e\ and e s . 

To solve Maxwell equations for the electromagnetic 
field operators in the Heisenberg's picture we introduce, 
in the usual manner (l3| . the electromagnetic potentials 
A(r,i) and $(r, i) and work in the generalized Coulomb 
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gauge 



V • [e(r)A(r)] = 0, 



(1) 



with the dielectric permittivity being a piecewise con- 
stant function as shown in Fig. [T] In the absence of free 
charges one can set $(r, t) — and work only with the 
vector potential A(r, t) which satisfies the wave equation 



V 2 A(r,i) 



d 2 



0, \z\?L/2. (2) 



Note that right on the interfaces condition (U} is singular 
due to discontinuities of the dielectric function and equa- 
tion ([2]) does not hold at these points. The normal-modes 
of the field f(r)e lult satisfy the Hclmholtz equation 



V 2 f kA (r) + e(z)w 2 f kA (r) = 0, \z\ ^ L/2, 



(3) 



and we have labelled them by their wave-vector k and 
polarization A = {TE, TM}. This mode decomposition 
allows one to solve the field equation (|5J) in each distinct 
region of space separately and then stitch up the solutions 
across the interfaces by demanding that they are consis- 
tent with the Maxwell boundary conditions, i.e. that E||, 
D±, and B are all continuous. 

The Hclmholtz equation (j3|) is in fact the eigenvalue 
problem of an Hcrmitean operator [12j 



V x V x 



VrfkA(r) = -w 2 yrf kA (r), (4) 



so that we expect the field modes \/efk.x(r) to form a 
complete set of functions suitable for describing any field 
configuration. The completeness relation takes the form 

|d 2 k||^ A (r)/*i(r') = ^.(r,r'), z,z' > L/2 (5) 



with &y(r, r') being the unit kernel in the subspace of 
functions satisfying ((TJ) ; we shall call this the generalized 
transverse delta-function. From quite general considera- 
tions [14| we can expect it to be given by 

^•(r, r') = «%«5 (3) (r - r') - G(r, r') (6) 

with the electrostatic Green's function of the Laplace 
equation given by 



G(r.r') 



1 1 



1 



47T |r - r'| 47r J 



dfcJ (M e~ k ( z+z "> 



n l ^ _ n l n s c -2kL 

n? + 1 n 2 + n 2 



nf-lnf-n 2 c _ 2kL 
n? + 1 n 2 + n 2 



(7) 
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where p = |r|| — r|i| and, for brevity, we have chosen to 
confine ourselves to the case z,z' > L/2. The function Jo 
in the above equation is a Bessel function of the first kind 



[Ta , 9.1.1]. The outline of the derivation of the Green's 
function is given in Appendix [B] 

The sum over all modes in equation §5§ is complicated 
because the spectrum of the field modes has non-trivial 
structure. It has been shown previously [l6| that 
the system supports two kinds of quite distinct types of 
modes. There are travelling modes going from left to 
right or in the opposite direction, and there are guided 
modes that are trapped by the dielectric layer, which 
essentially acts as a wave-guide. The spectrum of the 
travelling modes is continuous whereas the spectrum of 
the modes trapped in the dielectric layer is discrete and 
only some values of the (perpendicular) wave vector are 
allowed, namely those satisfying a certain dispersion re- 
lation. This dual character of the spectrum of the field 
modes is a major obstacle in working with these modes 
and calculating, e.g. the energy shift of an atom nearby, 
but an elegant solution to this problem has been devel- 
oped in [171 ] , whose basic idea we follow here. 

We choose the normalization of the mode functions 
v / efkA(r) according to the convention 



|d 3 r £ (z)f^(r).f k , v (r) 
5 xy d^(k-k') 



travelling modes 
8\\>S( 2 '(k\\ — k(,)(5fc z fe^ trapped modes 



(8) 



Then, the electric field E(r) = — dtA(r) expanded in 
terms of the normal-modes can be written as 

E(r) = iJ2 J d2k \\^2 V^ akAfkA(r)e '" kt +H '°' 

A k z 

(9) 

where H.C. stands for Hcrmitean conjugate. The photon 
creation and annihilation operators, aj^ and a k A, satisfy 
bosonic commutation relation 



. r^(k-k') 

AV | £( 2 )(k|| - k'«)6, 



(10) 



where the top and bottom of the RHS corresponds to the 
travelling and trapped photons, respectively. In order to 
be able to write out the electromagnetic field operators 
explicitly one needs to solve the eigenvalue problem (|3]) 
and determine the spatial dependence of functions f k A (r) 
so we turn our attention to this now. 



A. Travelling modes 

Before we work out the travelling modes, for further 
convenience, we introduce Fresnel coefficients for a single 
interface. For that we assume that a plane wave is travel- 
ling from a medium with refractive index rib to a medium 
with the refractive index n a , and that the interface is the 
z = plane. Then, the standard Fresnel reflection and 
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transmission coefficients are given by [lc 



ik? 



' TE 



'TM 



"-zb "'za 
^zb ~t~ ^za 
feb/«b — fea/^a ,ba 



iba 

I TE — 



2A: 



zb 



^zb + fcz 



(11) 



2fc zb /n a n b 



^zb/«b + ta/^a ' ™ 



j, T "/,;./ " ; , fczb/"-b + ^za/n-a 



where /c z ; are the components of the wave vectors per- 
pendicular to the interface in the medium i = {a, b}. 

The geometry of the problem (cf. Fig. [T]) naturally di- 
vides the space into three distinct regions. Consequently 
there are three wave vectors to be distinguished. The 
wave vector in vacuum (z > L/2) 

k ± = (fc x ,fc y ,±fc z ) = (k||,±fc z ), (12) 

the wave vector in the dielectric layer (\z\ < L/2) 

kf = (k x ,k y ,±k zl ) = (k||,±fc zl ), (13) 

and the wavevector in the substrate (z < —L/2) 

kf = (k X ,ky,±k ZS ) = (k||,±fc ZS ). (14) 



i.e. 



operators act on the factors of the type e 
eA(kj t ) = eA(V)e lk i r . Polarization vectors defined in 
such a way are normalized to unity provided all three 
components of the wave vector are real. However, they 
are not of unit length in the case of evanescent waves 
which have wave vectors with pure imaginary compo- 
nents. The spatial dependence of the mode functions is 
worked out requiring that each mode consists of the in- 
coming, reflected and transmitted parts that are joined 
together by standard boundary conditions across the in- 
terfaces, i.e. that E||, D± and B arc continuous. From 
this it is straightforward to derive that the travelling 
modes of the system incident from the left, normalized 
according to ([5]) , are given by 



eA(V) 
(27r)2n s 



ik+- 



1 A e 



z < -L/2 
\z\ < L/2 
z > L/2 



(20) 



whereas the right-incident modes are given by 



The components of the wave vector that are parallel 
to the surface are the same for all three regions of 
space. This follows directly from the requirement that 
the boundary conditions must be satisfied at all points 
of a given surface i.e. the spatial phase factors e lki ' r must 
be equal at z = ±L/2 for all rii. The different signs of 
the ^-components of the wave vectors correspond to the 
waves propagating in different directions. However, the 
direction of the propagation of a particular mode needs 
to be consistent in all three layers so we require that on 
the real axis 



sign(k z ) = sign(k z i) = sign(k zs ). 



(15) 



Since the frequency w of a single mode is fixed, the z- 
components of the wave vectors in the dielectric are re- 
lated to the vacuum wave vector k z by 



(n 



n\kl , 



fes = J{nl-l)k\+nlkl 



(16) 
(17) 



The mode functions ftA (r) are transverse everywhere ex- 
cept right on the interfaces z = ±L/2, cf. ((T|). To ensure 
this transversality, it is convenient to introduce orthonor- 
mal polarisation vectors 



f kA (r) = e A (k)/ kA (r) 



(18) 



defined as 



e TE (V) = (-A||)- 1 / 2 (-iV ys iV x ,0), 

e TM (V) = (A||A)- 1 / 2 (-V x V z ,-V y V z ,A||), (19) 



rpR ik s r 



(27r) 2 j , k .,. L jjH e *k + -r 



z < -L/2 
\z\ < L/2 
z > L/2 



(21) 



For the sake of clarity the complete list of reflection and 
transmission coefficients is given in Appendix [X] Here we 
only write down the ones most relevant for the calculation 
at hand: 



iff = 



rl 1 + r^e 



\s p 2ik z \L 



ik z L 



j-sl j-lv „(2ife z i -ifc ZB -ifc z )L/2 
rpL _ L \ L \ c 

* \ -\- ^.slj,lVg2ife zl L 

A A 



(22) 
(23) 



B. Trapped modes 

Trapped modes arise from repeated total internal re- 
flections within the layer of higher refractive index n\. 
This happens when the angle of incidence of the incom- 
ing wave is sufficiently high and exceeds the critical an- 
gle. This critical angle is different for the two opposite 
waveguide interfaces. First consider the layer- vacuum in- 
terface. From equation (|16p we can obtain the reciprocal 
relation expressing the k z in terms of the k z \ 



(24) 



Thus, whenever < (n 2 — l)k2 then k z becomes pure 



imaginary 



with A being the Laplace operator expressed in Carte- 
sian coordinates and it is understood that the above 



(n 



k 2 



(25) 
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and we have a mode that exhibits evanescent behaviour 
on the vacuum side. The sign of the square root is chosen 
such that these modes decay exponentially when one goes 
away from the layer in the positive z-direction. This also 
ensures that there truly is total internal reflection, i.e. 



that 



vl|2 



= 1. 



However, since on the other side of the waveguide we 
have a substrate rather than vacuum, not all of the modes 
that get totally internally reflected at the vacuum-layer 
interface necessarily get trapped. From the relation 



relation in terms of k z rather than k z \. It is straightfor- 
ward to show that the allowed values of the z-component 
of the evanescent waves' wave vector appearing on the 
vacuum side are given by numbers 



Ite = : fa + ik z \(k z ) tan[</>TE(£;z)] = °} - 
Qtm = {kz '■ K + z/c z i(fc z )/nf tan[</>TM(fc z )] = 0} 



with 



(33) 



k 

^ZS 1 

m 



(26) 



0te(A: z ) = arg [(k zl + k zs )e ifezlL ] , 
(f>TE(h) = arg [(fe z i/nf + k zs /n 2 s )e~ lk ° lL ] 



we obtain the condition of total internal reflection for 

.2 

'li- 



the substrate-layer interface to be k\ < (nf/n 2 — l)kjj 
Therefore, modes satisfying the condition 



(n 2 /n 2 



l)kf < kl < (n 



l)k 



(27) 



The numbers q x lie on the imaginary /c z -axis; they satisfy, 
cf. Eq. and (J28]), 



(34) 



-l-l)kf<(^) 2 <fi-l)ki 



ii: 



are not trapped but appear in vacuum as a continu- 
ous spectrum of evanescent waves that are accounted 
for among the left-incident travelling modes. (They are 
analogous to the evanescent modes that occur at a single- 
interface half-space, for which the normal-mode quanti- 
zation was first presented in Q.) On the other hand, 
trapped modes occur if 



< kl < (n 2 /n 2 



(28) 



The procedure for obtaining the trapped modes is largely 
equivalent to that of the travelling modes. They can be 
written in the form 



T]fe ik = - r 



z < -L/2 

f k ' A |ri = N x e x (V) { V x ^ r + e ik i + ' r \z\ < L/2 

z > L/2 



t^Iv 'ik -r 
A e 



(29) 



The boundary conditions are imposed on both interfaces. 
From the boundary at z = —L/2 we get 



jrls _ fj\si\s^ e -i(k^ l + k^)L/2 

whereas from the z = L/2 boundary 



/Tnlv _ Av„-i(k m \-kz)L/2 

2 x — l x e 



Vx 



— i \ c 



(30) 



(31) 



Since both equations, (|3"0")) and (|3II) . need to be simulta- 
neously satisfied we obtain a dispersion relation for these 
modes, 



1 + r v x l r^e 2iKiL 



0, 



(32) 



The normalization constant N\ for trapped modes is eas- 
ily obtained by direct evaluation of the integral © . It is 
given by 



Nx = ±- 
2ir 



1-1/2 



2n 2 L + F x (n h n s ) + F x (n h l) 



(35) 



with 



F x (n h n s 



|2 I^I 2 



--ilm(r^)e^(k+ 

K zl 



•eA(kj ) 



and the reader is reminded that in (1351) the z-components 
of the wave vectors k and k s are pure imaginary and 
because of that the TM polarization vectors eTivi(k~) 
and exivi(k~) are no longer normalized to unity, i.e. 
|e TM (k-)| 2 £ 1. 



C. Field operators and commutation relations. 
Completeness of the modes. 

Now that we have determined the spatial dependence 
of the mode functions we are in position to write out the 
vector potential field operator explicitly 

A(r,t) = { /d 2 k,| /"dfc.-^^f&Wafte-*** 



d'kn / dfe 



o V2eo^k 
I 



I d 2 kn Y, -=^C A (r)4 A e-^ k 4 + H.C. (36) 
J ^ V2e Wk J 



which determines the allowed values of k z \ within the 
layer. Since we will be dealing with an atom on the vac- 
uum side it will be necessary to express the dispersion 



The sum in the last term runs over the allowed values 
of the z-component of the layer's wave vector k z \, i.e. 
the solutions of the dispersion relation (|32p . For a given 
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type of mode, left-incident, right-incident, or trapped, 
photon creation and annihilation operators appearing in 
([36|) satisfy the commutation relations (fTTj|) . Commuta- 
tors between photon operators corresponding to different 
types of modes vanish as a consequence of the orthogo- 
nality of the field modes (jHJl, e.g. 



0. 



(37) 



We would like to verify explicitly the equal-time canon- 
ical commutation relation between field operators, say, 
between the electric field operator E(r, t) and the vector 
potential operator A(r, t) 



A l {r,t),e E J (r',t) = -iSUr, r'), z, z' > Z/2 (38) 



with 5y(r,r') given by Eq. ^ and (JTJ). To evaluate 
we shall need the electric field operator which is easily 
obtained from Eq. (|3f)t using the relation E = —d t A. 
Plugging in the field operators into (|3"5)) and making use 
of commutation relations (flU)) and (|3"T)l . we find that the 
LHS of (05]) is given by 



LHS = iReJ2 /d 2 k,| J^dk z f« Kl (r)fZl ;(*') 



+ I dk zs / k L A .,(r)/^,,(r') 



E 



/ k T v(r)/^(r') 



(39) 



The quantity on the right-hand side is the sum over all 
modes, just as prescribed by equation (O, and therefore 
we expect it to be equal to the generalized transverse 
delta function, Eq. ©. This shows that the statement of 
the completeness of the modes (0 is in fact equivalent to 
the commutation relation (|3"51) . as has been noted before 
in [l8| . To prove that the relation 

//■oo 
d 2 k|| / dk z f&/T)f£ j (T') 

p OO 

+ J dW k \,:(r)/k^(r') 



(40) 



holds for z, z' > L/2 we need to work out the sum over 
all field modes. To start with we carry out a change 
of variables in (j40|) : we convert the fc zs -integral and the 
fc z i-sum to run over the values of k z . In the case of the 
fc zs -integral this is a simple change of variables according 
to ([TiJ) 

/ dfc zs = n s 2 / dk z -^+n 2 s dfc z -^ (41) 

with T s = ^J(n 2 — l)k|/n s . Here it is seen explicitly that 
the contributions from the left-incident modes split into 



a travelling part and an evanescent part. The values of 
k z included in the last integral correspond to the con- 
dition for evanescent modes with continuous spectrum, 
Eq. (j2"7| . In the case of the sum we change the sum- 
mation over k z \ to run over the values of k z as defined 
by equation (|3"3")l . Plugging in the mode functions ([20)1 
and (j2"Tj) into equation (|40p and utilizing straightforward 
properties of the reflection and transmission coefficients 



that hold for real k z ,k zs , 



R\ R (-k z 



fez \r r L\ 2 



A 



(42) 



we can rewrite the completeness relation as 
^(r,r')=^(r-r') 



E^ 

X 



lj2\ N ^\ 2 K\ 2 elUz+z,) 



dky ^\Tt\ 2 e lk " (z+z,) 



(2tt) s 



dk z R^e lk ^ z+z "> 



(43) 



The first term in the above equation is the standard 
transverse delta-function. Therefore, if equation (|4U|) is 
to hold, the term in the curly brackets needs to be propor- 
tional to the reflection part of the electrostatic Green's 
function, cf. the second term on the RHS of Eq. 
That this is indeed the case is at this stage far from ob- 
vious, as for the proof one would need to combine two 
integrals and a sum into one expression. Obviously, the 
discreteness of the spectrum of the trapped modes is a 
nuisance that needs to be overcome if one is to complete 
the task of summing over the electromagnetic modes suc- 
cessfully. A similar difficulty would arise in any perturba- 
tive calculation in this type of geometry, which motivated 
a previous investigation of this problem for the symmet- 
ric case of a single slab of dielectric material ll7J. We 
proceed with a broadly analogous method to [17[, first 
noting that what we have here can be considered as a su- 
perposition of a slab and a half-space geometry, cf. [I?} 
and [IH . One can utilize the branch-cut due to k zs (which 
runs along the imaginary k z axis between ±ir s , cf. Fig. [2]) 
to express the integral over \T^\ 2 in (|4"3")l as an integral 
over the reflection coefficient that runs from 0~ along 
the square root cut up to the branch-point at +iT s and 
then back down to the origin + . Note that the branch- 
cut due to the k z \ is irrelevant because of the symmetry 
property of the reflection coefficient R^(—k z \) = R^(k z \). 
In this way, the first two integrals in the curly braces in 
equation (|4"3")) can be combined together as a single inte- 
gral in the complex k z plane [H[ . This is possible because 
the relation 



k ZB ,kzi>0 



-m 



fc ZB ,fc z i>0 



(44) 



fe zs ,fc a i<0 
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continues to hold for coefficients ([22]) with a purely imag- 
inary z-component of the vacuum wave vector, k z (cf. 
[19]). Thus, the contributions from the travelling and 
evanescent modes can be combined into a single con- 
tour integral along the path 7 S depicted in Fig. [2] and 
the terms appearing in the curly brackets on the RHS of 
Eq. become 



(27T) 



^3 J dfc z i?M(k+)ei(k-> 



ifc z (z+z') 



+ J2\Nx\ 2 K\ 2 i\(k + )e{(k-y k 



;(*+«') 



(45) 



Here we have now included the polarization vectors ex- 
plicitly in the integrals, which is a crucial step as they af- 
fect the analytical structure of the integrand in the com- 
plex A^-plane. In particular, the TM polarization vector 
introduces a pole at the points k z = ±£|kn | due to the fac- 
tor l/|k| 2 in its normalization factor. We will see that it 
is precisely this pole that gives rise to the reflection term 
in ©. We note that, according to Eq. (|2"21 . the reflec- 



Im(fc z ) 



Re(fc z 



FIG. 2: The dashed line represents the contour 7 S used to eval- 
uate the k z integral in Eq. (|45[) . Here r s = ^(n| — l)k|/n s 

and Ti = ^J(nf — l)k|/ni. The crosses represent the poles of 

the reflection coefficient i.e. the solutions to the dispersion 
relation (l32l). 



tion coefficient contains the phase factor 



-ik z L 



Thus, 



since z + z' — L > 0, the argument of the exponential in 
(|4"5| has a negative real part in the upper half of the 
complex k z plane and we can evaluate the /c z -integral in 
Eq. (|45|) by closing the contour in the upper half-plane. 
For this we need to determine the analytical properties 
of . We note that the denominator of the reflection 
coefficient (|2"2"j) is precisely the dispersion relation 



Rewriting the reflection coefficients in the form 



ky] 



-fi TE — 



1 — r!^ E exp(2ifc z iL) 
1 + rJ^ E exp(2ifc zl L) 



R R 



TM 



k 

A, z „ 

Tlf 



1 - rf E exp(2ifc zl L) 



1 + r}f E exp(2ifc z iL 

.Is 

TM 



1 - cxp(2ifc zl L) 



1 



'TM 



exp(2«fc z iL) 



fczl 



1 ~ r TM exp(2ifc z iL) 
1 + i 



' TM 



exp(2«fc z iL) 



allows us to deduce that Rf has a finite number of simple 
poles on the imaginary axis. When closing the contour 
we enclose all of them and by Cauchy's theorem the prob- 
lem is reduced to the evaluation of the residues at these 
points: 



£ / dk z R*e\(k+)e{(k-y k ^ + ^ 
x J f° 

= 2mJ2J2 *M ( k+ ( k ~ Y K (Z+Z ' } 

X Res 

V V lim (fc z - ql) + lim (fc z - i\ku 

* ' fe z — fr_— v«lk„l 



27TZ 



A 



; e\(k+)ei(k-) 



k z ~ >i|kii 
^.vl _|_ r ls e 2i/c zl L 
\ _|_ r vl r kg2ife z iL ' 



(46) 



Here, the first term represents the contributions from 
the poles in the reflection coefficient and corresponds the 
trapped modes, whereas the second term represents the 
contribution from the pole that arises due to the TM 
polarization vector. When calculating the residues ex- 
plicitly one needs to remember that the two independent 
variables are k z and kn and that, according to Eq. (fT6|) 
and (fT"?| . k z \ and k zs are functions of those. In addition, 
the denominator of the reflection coefficient is not of the 
form f(k z )(k z — g") so that multiplying it by (fc z — g") 
does not remove its singularity; the whole expression is 
still indeterminate. Therefore, L'Hospital's rule needs to 
be used to evaluate the limit (cf.[l7], Section V]). Doing 
so, we find that 

^3 J dk z R^i\(k + )i{(k-)e^ + ^ 

= - E l^l 2 \T{ v \ 2 e\(k + )e{(k-y k ^ + ^ 



-V 4 V'G ff (r,r') 



(47) 



where G#(r, r') is the reflected part of the Green's func- 
tion of the Poisson equation given in Eq. (|7|) and derived 
in Appendix [B] We see that the poles of the reflection 
coefficient i?f yield a term that exactly cancels out the 
contributions of the trapped modes to the completeness 
relation (14"5|) whereas the pole of the TM polarization 
vector yields the term proportional to Green's function. 



8 



Thus, the final result can be written as 



d 2 k|, TjL(r)&) = 7 [^(r),-£o^(r')] 



= <^(r - r') - V l V J G H (r, r') z, z' > L/2 
= %(5 (3) (r - r') - V,V;G(r, r') z, z' > L/2 

which is precisely what we have anticipated earlier. In 
the next section we demonstrate how the calculation pre- 
sented here may be applied to accomplish typical pertur- 
bative QED calculations in a layered geometry. 



III. ENERGY SHIFT 

To work out the energy shift we use standard pertur- 
bation theory where the atom is treated by means of 
the Schrodinger quantum mechanics and only the elec- 
tromagnetic field is second-quantized. We work with a 
multipolar coupling where the lowest order of the inter- 
action Hamiltonian is 



(48) 



Then the energy shift of the atomic state i, up to the 
second-order, is given by 



AE l = (i;0\H int \i;0)+Y, 




IC?;0|g in t|i;lkA)| : 
Ei - (Ej + w k ) 



Here, \i is the atomic electric dipole moment, and the 
composite state \j; describes the atom in the state 
\j) with energy Ej and the photon field containing one 
photon with momentum k and polarization A. Because 
the electric field operator is linear in the photon creation 
and annihilation operators, the first-order contribution 
vanishes and the second-order correction is the lowest- 
order contribution. Since the electric field does not vary 
appreciably over the size of the atom we use the elec- 
tric dipole approximation. Then the energy shift can be 
expressed as 



AEi = 




E i 



(49) 



where ro = (0, 0, zo) is the position of the atom and we 
have abbreviated Eji = Ej — Ei. It is seen that the 
calculation involves a summation over the modes of the 
electromagnetic field as carried out in the proof of the 
completeness relation (|4"3")l . Equation (|49p can be written 
out explicitly as 



2 /dk|| 

x (A vac + A trav + A cvan + A trap ) (50) 



with |/z m | 2 = |(*|Mm|j)| 2 - There are four distinct con- 
tributions to the energy shift. A vac is the position- 
independent contribution caused by the vacuum fields 
and gives rise to the Lamb shift in free space 



^ vac 



1 



(2tt) 3 



dk z e?(k-)er(k-) 



E 



(51) 
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The remaining three contributions come from the travel- 
ling, evanescent, and trapped modes, respectively, 



1 r°° 

AtraV = 7^3 / dfc ZJ R«er(k + )er(k-)e 



2ik z zo 



Eji + 



^evan 



(2tt) 3 



dfcz -^|T A L | 2 e^(k+)e^*(k+)e 2 ^ Z0 



Eji + LO 



(52) 



A trap = y-|^ A |2| r ^v|2 e m( k +) e m*( k +) e 2ifc zZ0 



Ei 



- OJ 



with Zq being the position of the atom with respect to 
the origin. Note that because of the dipole approxima- 
tion the shorthand notation for polarisation vectors (fT9|) 
can be no longer applied. Normally one is interested in 
the energy shift caused by the presence of the dielectric 
boundaries only i.e. the correction to the shift that would 
appear in the free space. Therefore, we renormalize the 
energy-level shift (|5U|) by subtracting from it its free space 
limit, i.e. 



AE, 



lim 

n\,n s — > 



AE, 



(53) 



The renormalization procedure amounts to the removal 
of the contributions A vac , Eq. ff5T|) - from the energy shift 
(|50p and takes care of any infinities that would appear 
otherwise, provided we treat the remaining parts with 
care. As noted elsewhere the contributions (|52l) suf- 
fer from convergence problems when treated separately. 
However, appropriate tools to handle the problem have 
been developed in Sec. Ill CI We aim to combine A trav , 
A ovan and A trap into one compact expression that is easy 
to handle analytically. We can use the same trick as in 
the proof of the completeness relation because the ana- 
lytical structure of the integrand in the complex A: z -plane 
is the same except for the function uj = (kjj + & 2 ) 1 / 2 
that comes about due to the denominator of perturba- 
tion theory and introduces additional branch-points at 
k z = ±i|kii | as compared to Fig. [5] This poses no difficul- 
ties though, if one chooses the branch-cuts to lie between 
±i|k||| and ±zoo. Then, the contributions to the energy 
shift from the travelling modes A trav and the evanescent 
modes A evan can be combined together into a single com- 
plex integral as explained in the steps between Eq. (|43|) 
and Eq. ([4"Sf . This is possible because for imaginary k z 
we have e™*(k + ) = e™(k~), whereas for real k z the re- 



lation e^*(k") 



l (k ) holds. On the other hand, we 



also know from Eq. (|4"T)) that the sum in A trap is equal to 
an integral over the reflection coefficient taken along 
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any clockwise contour enclosing all of it's poles. Choosing 
this contour to run from k z = 0~ + iT s to fc z = 0~ + iT\ 
and then back down from k z = + + iT\ to fc z = + + iT s , 
cf. Fig. [31 we write down the renormalized energy shift 
compactly as 



AE'™ = 



1 



2(2^) 3 e 



EE 



dk 



dfc z 



Co' 



<er(k+)er(k- 



(54) 



where the contour of integration 71 is shown in Fig. [3] It 
resembles that of Fig. [5] but now runs on the imaginary 
axis up to the point k z = iT\ enclosing all the poles of 
the reflection coefficients R^. Formula (f5"4"f is equally 



.'- 

1 

1 ) 


Im(fc z ) 
■ i|k||| 




i 

I 
ii 

i 

_ _J 


\iT s 

V. 






w- 

Re(fcz) 



FIG. 3: The dashed line represents the final contour 71 used 
to evaluate the energy shift in Eq. (154[) . 

applicable to ground-state atoms |0) as it is to atoms that 
are in an excited state \i) provided we use the contour 
of integration as given in Fig. [3] and interpret the k z 
integral as a Cauchy principal-value. As renormalization 
has now been dealt with we shall from now on omit the 
superscript "ren" and designate the renormalized energy 
shift of Eq. (JSIJ) simply by AEi. 



A. Ground state atoms 

In the case of a ground-state atom the energy differ- 
ence EjQ = Ej — Eq is always positive hence the de- 
nominator in Eq. (|54[) that originates from second-order 
perturbation theory, Ejq + uj, never vanishes. Then, Eq. 
([54)) contains no poles in the upper half of the fc z -plane 



other than those due to the reflection coefficient R x . To 
evaluate the k z integral we can deform the contour of 
integration in Eq. ((54)) from that sketched in Fig. [3J 
to the one as shown in Fig. @] which is beneficial from 
the computational point of view as it simplifies the anal- 
ysis of Eq. (fMf considerably. Writing out explicitly 



!lm(fc z ) 
1 

Hk||l 



Re(k z ] 



FIG. 4: The final contour C used to evaluate the energy shift 
of the ground state atom in Eq. (|55[) . 



the sums over the polarization vectors (119)) and then ex- 
pressing the integral in the ky -plane in polar coordinates, 
fc x = fcii cos <f>, k y — fcii sin <fi, where the angle integral is 
computable analytically, we rewrite the energy shift as 

1 v 



167r 2 eo , 



BR 



TE 



dfc|| 



k 2 ~ 

z £> R 
~ 2 n TM 



dfc z 



„2ik z Z 



' E 



jo 



-2|mj 



I pi? 



TM 



(55) 



with u(k z ) = ^//c^ + /cj, \fi\\\ 

contour C is that in Fig. 2) 
coefficients R^ are given by 



2 and the 



= IMzI + 

The amended reflection 



Is 2ik zl L 

r x e 



1 + r v ^r^ s e 2ikilL ' 



(56) 



i.e. we have pulled out the phase factor e~ lklL in order 
to define Z = zq — L/2 as the distance between the atom 
and the surface, cf. Eq. (12^1) . 

In order to perform the fc z integration in (|55|) we need 
to analytically continue the function uj = w(fc z ), which is 
real and positive on the real axis, to the both sides of the 
branch cut along which the integration is carried out, (cf. 
Fig. [4j. Doing so we find that on the LHS of the cut the 
positive value of the square root needs to be taken, and 
hence on the RHS of the cut we must take the opposite 
sign. Therefore we have 



dfc z 



dfcz 



lEjQLJ 



1Q 



(E j0 - u)(E j0 +u>) 



Now we carry out a sequence of changes of variables. 
First we re-express the k z integration in terms of one 

over the frequency uj by substituting uj - 



fc 2 1 



dfcz = 



duj 



uj 2 - fcjj 



(57) 
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Then, we make the integral run along the real axis by 
setting uj = i£. After this is done, the energy shift of the 
ground state is expressed as a double integral that covers 
the first quadrant of the (k\\ , £)-plane 



AE = 



{hl 2 [« 2 



Ejo I dk\\ fc|| 



kf\ )^TM 



d£ 



C 2 f>R 



It seems natural to introduce polar coordinates, fen = 
xsin0, £ = 5cos0. We also choose to scale the radial 
integration variable S = -Ejo^ with Ejo > and set y = 
cos 0. This provides us with the final form of the energy 
shift that is more suitable for numerical computations 
and asymptotic analysis 



AE 



— Y 



E 



jo 



dxx c 



dy 



-2EjqZx 
2, ,2 



o jo l + x-y 

Jt-' 

lMiil 2 (y 2 ^TE-^TM, +2\v±\ 2 (y 2 - ll/? 



R 

TM 



(58) 



The reflection coefficients are as expressed in ([56)) but 
with the wave vectors given by 



fc zi = ixE j0 yJ (n 2 - l)y 2 + 1, Th = {1, ni, n s }. 

Note that even though the wave vector is imaginary, the 
final result is a real number, as it should, because the 
Fresncl coefficients contain only ratios of wave vectors. 



B. Excited atoms 

As mentioned previously, the energy-level shift of an 
excited atom is also given by Eq. (|54p. However, one 
needs to take account of the fact that the quantity 
Eji = Ej — Ei can now become negative for Ej < E.i, 
so that the denominator originating from perturbation 
theory contributes additional poles lying on the path of 
k z integration, shown in Fig. [3] and is now to be under- 
stood as a Cauchy principal-value. These poles are lo- 
cated at k z = ± yJlZji — kjj , though their precise location 

depends on the value of |k|| | that is not fixed but varies 
as we carry out the k|| integrations in equation (|54|) . For 
|k||| € [0, \Eji\] the poles are located on the real k z axis 
but as we increase the value of |k|| | to exceed \Eji\ both 
poles move onto the positive imaginary axis according to 
the convention that Im(k z ) > 0. For |k|| | belonging to the 
interval [|-Ej'i|, risl-E^'i |] the poles are located on the oppo- 
site sides of the branch-cut due to the k zs and care needs 
to be taken when evaluating those pole contributions. To 
evaluate the Cauchy principal-value of the /c 2 -integral we 
circumvent the poles and close the contour in the upper 
half-plane, as was done in the previous section. The con- 
tribution from the large semicircle vanishes and equation 



(|54p acquires pole contributions that are easily worked 
out by the residue theorem. The energy shift splits into 
the a "non-resonant" ground-state-like part AEi and a 
"resonant" oscillatory part A£"[ es that arises only if the 
atom is in an excited state. In analogy to the result of 
the previous section, the "non-resonant" part is given by 



e + k 2 (E% + 

+ 2k\R% M \ii 1 _ 



AEi = 



1 



i 3 } 



. rl e -2\ Eji \Zx 

dxx dy — , 
Jo l + x 2 y 



2, .2 



M||| 2 {V 2 Rte -Rtm) + %±lV - 1)^ M J (59) 
with wave vectors expressed as 

fc zi = ix\Eji\^J (n 2 - l)y 2 + 1, rii = {1, n h n s }, (60) 
whereas the "resonant" part is given by 
AEl cs = Re— Y \E rt \ 3 f°° dqQ P 2 ^\V^z 

87T60 ^ J 



3<i 

„2\ f,R 



x{|M||| 2 [(l-q Z )R^ M -R TE 

with wave vectors expressed as 



2|^|V^™},(61) 



fczi = {Ejil^Jn 2 - q 2 , Ui = {l,n h n s }. 

The reflection coefficients are as given in (f56|). The 
integral in Eq. (|6"Tj) contains poles because the disper- 
sion relation present in the denominators of the reflec- 
tion coefficients has now solutions on the real axis when 
q € [n s ,ni]. This signals contributions from surface exci- 
tations (trapped modes). This fact has been mentioned 
in [2| where the interaction of an excited atom with lay- 
ered dielectric has been studied, although using mainly 
numerical analysis. Here we will attempt to study the 
results (|59p and (|6"Tj) analytically. To do so it will prove 
beneficial to rewrite equation (|6"Tj) slightly. We change 
variables according to \/l — q 2 = r\ and split the contri- 
butions to Eq. (|6ip into two parts. The first one is a 
contribution from the travelling modes and given by 



AE] 



-Re 



2i\Eji\Zr) 



x {hi 2 [r^ e -v 2 Rtm]+^±\ 2 (1-v 2 )Rtm} (62) 

where the wave vectors in reflection coefficients are all 
real and can be expressed as 



fczi = \Eji\ynf - 1 + if, rii = {l,n h n s }, (63) 

and the second is a contribution from the evanescent 
modes 



AE, 



res, e van 



-Re- 



1 r c 
— T \En\ 3 \ 



drje 



-2\Eji\Zri 



{hi 2 [r^ e + v 2 Rtm]+2\^\ 2 (1+V 2 )Rtm} (64) 
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where the wave vectors in reflection coefficients can be 
expressed as 



Ki = \Eji\sJnt - 1 - rf 



= {l,ni,n s }. 



(65) 



Finally, it is worth noting that the imaginary part of Eq. 
(l6"Tj) is actually proportional to the modified decay rates 
Q. These have already been studied in 16] so that we 
focus on energy shifts only. However, the methods of 
analysis that are reported in the next section do allow 
one to write down at once equivalent analytical formulae 
for the decay rates. 



IV. ASYMPTOTIC ANALYSIS 

The interaction between the atom and the dielectric is 
electromagnetic in nature and it is mediated by photons. 
The atomic system in state \i) evolves in time with a char- 
acteristic time-scale that is proportional to Ej t , with Eji 
being the energy- level spacing between the states \i) and 
|j) which are connected by the strongest dipole transi- 
tion from state \i). Since it takes a finite time for the 
photon to make a round trip between the atom and the 
surface, the atom will have changed by the time the pho- 
ton comes back. Therefore, the ratio of the time needed 
by the photon to travel to the surface and back and the 
typical time-scale of atomic evolution is a fundamental 
quantity that plays decisive role in characterising the in- 
teraction. In natural units, if 2EjiZ C 1 we can safely 
assume that the interaction is instantaneous and we are 
in the so-called non-retarded or van der Waals regime. If 
2EjiZ ^> 1 the interaction becomes manifestly retarded 
as the atom will have changed significantly by the time 
the photon comes back. However, the problem we have 
considered here provides us with yet another length scale, 
namely the thickness of the top layer L. We shall now 
consider the energy shift in various asymptotic regimes. 



A. Ground state atoms. Electrostatic limit, 

{2E ji Z « 1) 

In this limit the interaction is instantaneous (or elec- 
trostatic) in nature and the energy shift is obtainable 
using the Green's function of the classical Laplace equa- 
tion (cf. e.g. [20(). This classical derivation is outlined 
in the Appendix [Bj The end result for the energy shift 
reads 



AE cl = -- 



1 



167TCO 



{(4 



2 (iA 

n\ - 1 
n? + 1 



dfcr 



-2k Z 



's „-2fcL 



1 - 



In? 



-2kL 



,(66) 



with (/if) = (jil) + (fj, 2 y ) and (fi 2 ± ) = We will now 

show that one can also obtain the above result as a limit- 
ing case of the results of previous section, thus providing 
a cross-check for our general calculation. To start with we 
note that equation ([58| cannot be used to take the elec- 
trostatic limit in which we mathematically let Eji — > 
because it has been scaled with Eji. Therefore, it is best 
to start from equation (|54"|) . The result of Eq. (|6"6"]l can 
be derived very quickly if we observe that in the limit 
the branch cut due to u> = .. Ik 2 + k 2 is no longer 



present and the contour in Fig. 2] collapses to a simple 
enclosure of the point k z = i|kii|. The contribution from 
the TE mode vanishes as the product of the polarization 
vectors is regular at fc z = i|kii|, but for the TM mode this 
point is a simple pole, cf. Eq. ([!§]) . Therefore we obtain 



AE C[ = - 



1 



(27r)32e 



x27ri lim (k z 
k„-+i|k[i| 



EEi^i 2 / dk n 

m jjti 

i|k|||)^ M e^ M (k+)e^ M (k-)e 2 ^ 2 °. 



Taking the limit and expressing the remaining integrals in 
polar coordinates, where the angle integral is elementary, 
yields equation (g>6j with (tf n ) = \{i\fx m \j)\ 2 = 

(*l/" t mK)' Equation (|66p can be further analysed depend- 
ing on the relative values of L and Z. 



1. Thin layer (Z/L 2> I) 

In this case the distance of the atom from the surface is 
much greater than the thickness of the layer of refractive 
index n\ (but still small enough for the retardation to be 
neglected). Then, rescaling the integral in equation (166[) 
with k — x/L allows us to use Watson's lemma [22| to 
derive the following result 



AE C 



AEZ 



1 



647re Z 3 



((Mf> + 2(Mi> 



L L 2 fL 3 
ai -+a 2 — 2 +0\— z 



(67) 



with the coefficients aj given by 

o 4 4 

Oi 



«2 



n\ [nl + l) 2 ' 

6 (nf — nf)(n 2 + nf) 
K 2 + l) 3 



where AE° l is the well-known electrostatic interaction 
energy between an atom and a dielectric half-space of 
refractive index n s that can be obtained by the method 
of images 



AEl 



647T6 



(68) 



1 n{ 



The corrections to this result are represented by the re- 
maining elements of the asymptotic series. Note that if 
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n\ > n s then a\ > and, not surprisingly, the interac- 
tion, as compared to a half-space alone, is enhanced by 
the presence of the thin dielectric layer of higher refrac- 
tive index n\. 



distance of the atom from the surface and therefore is 
not measurable. The next-to-leading order correction is 
the first to be distance-dependent and is proportional to 
Z / L , which can be easily seen by expanding the factor 
in series around ZjvL = 0: 



2. Thick layer (Z/L^.1) 

In this case the thickness of the layer is much greater 
than the distance between the atom and the surface. The 
top layer now appears from the point of view of the atom 
almost as a half-space of refractive index n\ only that it 
is in fact of finite thickness. To analyse the result (j6"6")l in 
this limit we cast it in a somewhat different form. Note 
that, especially when kL is large but not only then, 



In 2 



-2kL 



In 2 



< 1 



(69) 



and the denominator of the integrand in Eq. (|66j) can 
be written as geometrical series. Since the series is abso- 
lutely convergent we can integrate it term by term and 
obtain the following representation of the electrostatic 
result 



AE C 



AEl 



1 



167rer 



<A*f> + 2</iI> 



E 



1 n\ 



In 



(Z + vLf 



(70) 



where AE^ is the electrostatic energy shift due to a sin- 



gle half-space of refractive index n\, i.e. Eq. (|68|) with n s 
replaced by n\. The sum in Eq. (|70p represents the cor- 
rection to AE~ due to the finite thickness of the layer. 
For fixed Z and L it can be easily computed numeri- 
cally to any desired degree of accuracy. We note how- 
ever, that to the leading order in Z/L the interaction 
is weakened by the same amount independently of the 



1 



1 



(Z + vLf v 3 L 3 



3Z 



(71) 



B. Ground state atoms. Retarded limit, 

(2ZEji > 1) 

1. Thin layer (Z/L > 1) 

In this case we study the situation when the top layer 
is much thinner than the distance between the atom and 
the surface. To obtain the asymptotic series we use Wat- 
son's lemma in much the same way as in the electrostatic 
case [2l|. Series expansion of the integrand in Eq. (|58|) 
about x — decouples the integrals and the resulting 
integrals can be calculated analytically. Thus, to first 
approximation, for an atom located sufficiently far from 
the interface, the impact of the thin dielectric layer on 
the standard Casimir-Polder interaction can be described 
by 



AE rct = AE i 
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(72) 



where A£"Jj et is the retarded limit of energy shift as 
caused by a single dielectric half-space of refractive in- 
dex n s , which was calculated in Q. We give this result 
in Appendix [C] The coefficients an and a± in (|72p can 
be expressed in terms of elementary functions as 
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Both, an and a±, are positive for n\ > n s so that, as 
one would expect, the interaction, as compared to a half- 
space alone, is enhanced by the thin dielectric layer of 
the higher refractive index n\. The above result simpli- 
fies significantly in the case when n s approaches unity i.e. 
when the situation resembles that of an atom interacting 
with a dielectric slab of refractive index n\. The coeffi- 
cients a || and a± reduce then to those recently calculated 
in [l[ and are given by 



l)(9n? + 5) 



10nf 
l)(5n 2 



4) 
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2. Thick layer (Z/L^l) 

Here we assume that the thickness of the top layer is 
much greater than the distance between the atom and the 
surface, but which is still large enough for retardation to 
occur. Note that the reflection coefficient R^ ([^| can 
be separated into L-dependent and L-independent parts 
in the following manner: 
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This way of writing the reflection coefficient splits the 
energy shift (f58|) into a shift due to the single interface 
of refractive index n\ and corrections due to the finite 
thickness and the underlying material. It can be shown 
numerically, see Sec. [V] that for large values of L the 
correction term is vanishingly small and can be safely 
discarded. Brute-force asymptotic analysis allows us to 
draw similar conclusions as in the electrostatic case, Sec- 
tion |IVA~2] To leading order the interaction gets altered 
by the same amount regardless of the position of the atom 
with respect to the interface. The next-to-leading-order 
correction is proportional to Z/L 5 . 



C. Excited atoms. Non-retarded limit, 

(2Z\En\ < 1) 



contributes only positive powers of Z. This is most eas- 
ily seen by expanding the exponential exp(2i\Eji\Zrf) 
about origin as we may do in the limit 2Z\Eji\ — > 0. 
Therefore, to leading-order in the electrostatic limit, only 
(f6"4"|) contributes. Further we analyse (|64p by setting 
■q = (3/(\Eji\Z). Then, according to (|65|) . in the limit 
\Eji\Z — > the wave vectors can effectively be approxi- 
mated as 



^ZS ~ ^ g 
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Then the result for the energy shift, after substituting 
f) = kZ, reduces to 
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This result turns out to have the same dependence on 
Z and L as the Coulomb interaction of the ground state 
atom, cf. Eq. (|6"6")l ; therefore we shall not analyse Eq. 
(f75|) any further. Note however, that the dependence 
on the atomic states is different in equations (|66[) and 
([75^1 . We would also like to point out that in the electro- 
static limit, to the order we are considering, the quantity 
AE res,el turns out to be real, which would imply that 
the corrections to the decay rates vanish. However, this 
conclusion is incorrect as it is known that the change 
of spontaneous emission in the non-retarded limit is in 
fact constant for a non-dispersive dielectric half-space |9j . 
However, any serious analysis of the changes of the de- 
cay rates induced by a surface needs to take into account 
the absorption of the material, which in the non-retarded 
limit plays a crucial role and cannot be neglected. Fur- 
thermore we note that we have started from Eq. (|61|). 
which, as explained before, contains poles on the real 
axis signalling the trapped modes. However, the denom- 
inator of (|75|) never vanishes which reflects the fact that 
in the electrostatic limit the trapped modes cease to exist 
and do not contribute towards the energy shifts, as first 
mentioned in 



The energy shift of an excited atom is given by equa- 
tions ([591 and (|6Tj) . The "non-resonant" part, i.e. Eq. 
(|S"5|) has the same form as the energy shift of the ground 
state atom and has been analysed in the previous sec- 
tion. Therefore we now focus on the "resonant" part of 
the interaction that is given by equation (1611) . In order to 
conveniently obtain the non-retarded limit of (|61|) we will 
work with its slightly modified form given in equations 
(HJ and (HII). 

We start by noting that close to the interface we expect 
asymptotic series to be in the inverse powers of Z. Equa- 
tion where the r\ integration runs over r\ £ [0,1], 



D. Excited atoms. Retarded limit, (2Z\Eji \ S> 1) 

The leading-order behaviour of equation (|61j) in the re- 
tarded limit can be obtained by repeated integration by 
parts. Unlike in the electrostatic case now both equa- 
tions, Eq. (|62|) and Eq. (|64|) contribute. We integrate 
them by parts and note that the non-oscillatory contri- 
butions that arise from the boundary terms evaluated at 
rj = cancel out. It turns out that the leading-order 
contributions to the energy shift are due to the perpen- 
dicular component of the atomic dipole moment. They 
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dominate the retarded interaction energy and behave as 
Z _1 . The contributions due to the component of the 
atomic dipole moment that is perpendicular to the sur- 
face contribute only terms proportional to Z~ 2 . We find 
that in the retarded limit the interaction energy up to 
the leading-order is given by 
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1 + 2r v m s cos(2\Eji\ T ) + r^r 2 
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+ri s cos[2\E ji \(Z + T)}}, (76) 

where we have defined the optical thickness of the layer 
as r = n\L and 



r v i 



i — m 
I 



ni - n s 
m + n s 



(77) 



The final result agrees with that derived for a half-space 
in Q if we take either L — > or n\ — > n s , which is a 
consistency check of our calculation. However, the limit 
of perfect reflectivity of the top layer does not make sense 
and one has to start from equation (|61[) and rewrite the 
reflection coefficient in the form (|75|) in order to study 
this case. 

Equation (|76[) is valid only approximately when the 
distance between the atom and the surface is much 
greater than the wavelength of the strongest atomic 
dipole transition, but it nevertheless allows us to draw 
important conclusions. We note that the interaction is 
resonant i.e. it is enhanced for certain values of LEji. 
The most convenient way to understand the essence of 
these resonance effects is to take the slab limit of equa- 
tion (|T6]) i.e. set n s — 1. In this limit we have 



AE* 
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complications arise from the fact that when the atom is 
close to the surface the evanescent waves come into play 
whereas the condition (f79| refers to the interaction of an 
atom with travelling modes only. In the non-retarded 
limit Z\Eji\ <C 1 the notion of resonance loses its mean- 
ing altogether, cf. Eq. (|75[) . Exploring the extreme case 
in the retarded limit we note that at anti-resonance i.e. 
when 



nL = —H-K, 
2 ' 



k = 0,1,2 .. 



(80) 



equation ([76]) becomes 
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V||| 2 cos(2|^|Z), 



(81) 

i.e. the atom does not feel the presence of the layer and 
the interaction assumes the form of that between an atom 
and a single half-space of refractive index n s , cf . Q . This 
means that in the retarded regime the leading-order in- 
teraction between an excited atom and a slab of thick- 
ness L vanishes whenever the optical thickness of the slab 
r = n\L is equal to a half-integer multiple of the wave- 
length of the dominant atomic transition Xji (cf. also Fig. 
QT] later on). Conversely, at resonance the shift becomes 
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V |^| 2 | M || I 2 cos^lZ), 



so that the amplitude of oscillations exceeds the ampli- 
tude that would have been caused by a single half-space 
of refractive index n\. It also reaches the perfect reflec- 
tor limit n\ — > oo more rapidly. Finally, we shall also 
remark that the meaning of the conditions (|79|) and (|80|) 
is interchanged if the refractive index of the substrate n s 
exceeds that of the layer n\ i.e. when n s > n\. 



V. NUMERICAL EXAMPLES 



l-2rJ 1 cos(2|£?ji|r) + < 1 
x {r vl (l + r^) cos(2|^|Z) 



cos[2|^|(Z-r)] 
-r vl cos[2|^|(Z + r)]}. (78) 



It is easily seen that whenever cos(2|i?jz|r) = 1 then 
AEl° s ' rc — 0, i.e. the leading-order interaction vanishes. 
Conversely, the amplitude of oscillations in equation (|78|) 
is maximized when cos(2\Eji\r) = — 1. Therefore we 
have a condition for resonance in terms of the wavelength 
of the strongest atomic dipole transition Xji 



nL 



X 



0,1,2 



(79) 



Eq. (J75J holds for Z\E^\ » 1 but if the value of Z\E^\ 
approaches unity, the relation loses its validity, because 



In this section we present a few numerical results de- 
signed to illustrate the influence of the dielectric layer on 
the Casimir-Polder interaction between an atom and a 
dielectric half-space. In practice, the sum over interme- 
diate states j in Eq. (|55|) and in Eq. (|ST1l is restricted to 
one or a few states to which there are strong dipole tran- 
sitions. Hence, we assume a two-level system in which 
Eji is a single number, namely the energy spacing of the 
levels with the strongest dipole transition. Additionally, 
we focus just on the contributions to the energy shift due 
to the component of the atomic dipole that is parallel to 
the interface of the dielectrics. The contributions due to 
the perpendicular components of the atomic dipole mo- 
ment can be easily generated with from Eq. (|58[) using 
standard computer algebra packages like Mathcmatica 
or Maple. We start by simple checks on the asymptotic 
expansions derived in the previous section. 
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A. Ground-state atoms 
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[5]and Fig. [6] We remark that even though the derivation 
of the energy shift in this paper was based on the assump- 
tion n\ > n s , the results are also valid in the case when 
the top layer has a smaller reflectivity than the substrate. 
In such a case the result can be used e.g. to model a thin 
layer of oxide or any kind of dirt on the substrate which 
is often present under realistic conditions. 

The asymptotic expansion ([7^]) works well for large 
Z/L and not too high values of the refractive index n\. 
This is demonstrated in Fig. [3 The increase of the 
refractive index n\ has an impact on the accuracy of the 
approximation which is valid provided 

Z > \ 3i + n (83) 

with Xji being the wavelength of the dominant atomic 
transition and r\ = n\L is the optical thickness of the 
top layer. In Fig. [8] we demonstrate the behaviour of 



FIG. 5: Plot of the exact energy-level shift contributions AE'' 
(solid), Eq. ([58j , multiplied by Z 4 . Dashed lines represent 
the energy shifts due to the single dielectric half-spaces of re- 
fractive indices n\ (top) and n B (bottom), whereas the dotted- 
dashed lines represents the asymptotic approximation (I72[) , 




FIG. 6: Plot of the exact energy-level shift AE" (solid), Eq. 
(|58|) . multiplied by Z 4 . Dashed lines represent the energy 
shifts due to the single dielectric half-spaces of refractive in- 
dices ni (bottom) and n a (top), whereas the dotted-dashed 
lines represents the asymptotic approximation (72) . 
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FIG. 7: Plot of the exact energy shift AE", (solid, Eq. ([58)} ), 
multiplied by Z 4 together with the asymptotic approxima- 
tions (dashed, Eq. ([72)) h 

the energy shift depending on the various values of the 
parameter Eji measured in units of the layer's thickness. 
For small Eji we clearly observe linear behaviour that 
corresponds to the Z~ 3 dependence of the shift in the 
electrostatic regime. 

We also find it instructive to plot the energy-level shift 
as a function of the thickness of the top layer L for dif- 
ferent values of the refractive index n\ while keeping the 
distance of the atom from the surface fixed, Fig. [9] and 
Fig. [HE 



We choose to plot the energy-level shift AE multiplied 
by Z so that the asymptotic behaviour of it as a func- 
tion of distance is more apparent, because Z 4 AE for a 
dielectric half-space approaches constant @. Then, one 
can easily track the variation of the energy shift caused 
by the top layer as compared to the half-space shifts, Fig. 



B. Excited atoms 

The energy shift of an excited atom splits into two 
distinct parts, cf. Eq. (I5T)1) and Eq. (|aTj) . The non- 
oscillatory part displays the same behaviour as the en- 
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FIG. 8: Plot of the exact energy shift AE^ (Eq. l|58)l) mul- 
tiplied by Z 4 as a function of Z/L for various values of the 
retardation parameter EjiL. 



FIG. 10: Plot of the exact energy shift AS 11 (Eq. lfg8]l) mul- 
tiplied by _> 4 as a function of layer's thickness L measured in 
units of fixed atom- wall separation Z for various values of the 
substrate's refractive index n s > n\. 
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FIG. 9: Plot of the exact energy shift AE n (Eq. J58J)) multi- 
plied by Z 4 as a function of layer's thickness L measured in 
units of fixed atom-wall separation Z for various values of the 
layer's refractive index n\ > n a . 



ergy shift of the ground-state atoms, which we have al- 
ready analysed numerically in the previous section. Here 
we will focus on the oscillatory contributions to the level 
shifts that are given by Eq. (l6"Tj) . We choose to plot the 
dimensionless integrals contained in equations (|62l) and 
(|64[) as this is numerically more efficient than plotting 
the integral in Eq. (|6ip . It should be borne in mind 
that the reflection coefficients contain the dispersion re- 
lation in denominators that now has solutions on the real 
axis. For the purpose of the present demonstration it is 
sufficient to simply displace the poles off the real axis by 
adding small imaginary part to the denominator of the re- 



flection coefficients, which amounts to taking the Cauchy 
principal-value during numerical integration. In Fig. 
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FIG. 11: Plot of the exact energy-level shift (|61[) (resonant 
part) in an excited atom due to the parallel component of the 
atomic dipole moment placed in front of a slab of thickness L 
and refractive index n\ = 2n. The energy spacing of the dom- 
inant atomic transition is such that LEji — 3/4 i.e. it satisfies 
the resonance condition (|80|) . As is seen, when LEji == 1/2, 
the energy shift in the retarded regime is strongly suppressed, 
cf. Eq. (EU). 

[TTlwe demonstrate that indeed, if the anti-resonance con- 
dition (|80)) is satisfied, the interaction energy between 
the excited atom and the slab is strongly suppressed for 
ZEji 3> 1. In general, for the layered dielectric rather 
than the slab, the effect of resonance is shown in Fig. _____ 
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FIG. 12: Plot of the exact energy-level shift (|61[) (resonant 
part) in an excited atom due to the parallel component of the 
atomic dipole moment placed in front of the layered dielectric 
with parameters as shown on the graph (solid). The resonant 
condition (|79[) is satisfied so that the interaction is enhanced. 
The amplitude of oscillations exceeds the one that would have 
been caused by an unlayered half-space of the refractive index 
n = 2-7T, cf. Eq. (|81[) . Compare also Fig. [5] The dashed lines 
represent the interaction between an atom and single half- 
space of refractive index n as indicated. 



and Fig. [T3l Note that the energy-level shift in an excited 
atom due to the layered dielectric can be significantly 
enhanced. Unlike in the case of the ground state atom 
where the energy shift caused by the layered structure 
of refractive indices n\ and n s is bounded by the single 
half-space shifts (compare Fig. [5]), the excited atom can 
experience shifts greater than those caused by the unlay- 
ered half-space of the refractive index n — max(m,n s ), 
Fig. [12j which is due to resonance effects. Conversely, 
it is also possible that the interaction with the layer will 
be unnoticeable if the anti-resonance condition (|80[) is 
satisfied, Fig. [13l Next, in Fig. [HI we show that the 
approximation of Eq. (|61[) derived in (|76[) turns out to 
be quite accurate and can be safely used to quickly esti- 
mate the energy shift in an excited atom caused by the 
layered dielectric, provided the condition ZEji ^> 1 is 
satisfied. It is also interesting to plot the resonant part 
of the energy shift as a function of LEji while keeping 
ZEji fixed. This is done in Fig. [TH It is seen that the 
energy shift indeed experiences the oscillatory resonant 
behaviour. The subsequent minima and maxima are less 
and less pronounced as the value of LEji increases. This 
is because as we increase LEji the resonances and anti- 
resonances move closer and closer together so that their 
effects cancel out. It is interesting to note that this be- 
haviour could not have been inferred from equation (|76p , 
which indicates that the approximation (|76|) can be use- 
ful only for LEji "C 1, which can also be easily verified 
numerically. 



FIG. 13: Plot of the exact energy-level shift (|61[) (resonant 
part) in an excited atom due to the parallel component of the 
atomic dipole moment placed in front of the layered dielectric 
with parameters as shown on the graph (solid). The anti- 
resonant condition (|80[) is satisfied so that the presence of 
the layer is almost unnoticeable, cf. Eq. (|81[) . The dashed 
lines represent the interaction between an atom and single 
half-space of refractive index n as indicated. 



VI. SUMMARY 

Using perturbation theory we have calculated the 
energy-level shift in a neutral atom placed in front of 
a layered dielectric half-space, as shown in Fig. [TJ The 
major difficulty in working out the energy shift is the 
sum over all modes that appears in this type of calcu- 
lation, Eq. (|50|) . especially when the spectrum of the 
modes consists of the continuous and discrete parts, Sec. 
UTAl and DIH This obstacle can be circumvented by us- 
ing complex- variable techniques to express the sum over 
all modes as a single contour integral in the complex k z - 
plane, Eq. ([54]) and Fig. [4] Then, the energy shift flSHJ is 
easily analyzed asymptotically as well as numerically. For 
a ground-state atom, regardless of whether in retarded or 
non-retarded regimes, we find that the leading-order cor- 
rection to the interaction of an atom with an unlayered 
interface is proportional to LjZ. The asymptotic series 
are given by (j67|) and (1721) and provide reasonable es- 
timate of the influence of the single dielectric layer on 
the standard half-space result, Fig. [7] In the opposite 
case of a very thick layer i.e. Z/L « 1 we find that 
the result is well approximated by a dielectric half-space 
Q. For excited atoms we find that the interaction be- 
tween an atom and the layered dielectric (j6Tj) is subject 
to resonances that occur between the wavelength of the 
dominant atomic transition Xji and the thickness of the 
layer L, Sec. IIVDI In particular, the interaction be- 
tween an atom and the slab can be strongly suppressed 
in the retarded regime, cf. Fig. 111! whenever the optical 
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FIG. 14: Plot of the exact energy-level shift (|61[) (resonant 
part) in an excited atom due to the parallel component of the 
atomic dipole moment placed in front of the layered dielectric 
with parameters as shown on the graph (solid). The dashed 
line represents the approximation in the retarded regime, Eq. 



tarded (|76|) regimes that can be used to quickly estimate 
the magnitude of the resonant interaction between an 
atom and a layered dielectric. 



Appendix A: Fresnel coefficients for layered 
dielectric 



Here we list the reflection and transmission coefficients 
appearing in the normal-modes of the system as discussed 
in the section Hi Al For the left-incident modes we find 
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and for the right-incident modes we get 
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The Fresnel reflection coefficients r^ b for a single interface 
are given by (fTTj) . 



Appendix B: Electrostatic calculation of the 
energy-level shift in a ground-state atom in a 
layered geometry 



FIG. 15: Plot of the exact energy-level shift (|61[) (resonant 
part) in an excited atom due to the parallel component of the 
atomic dipole moment placed in front of the layered dielectric 
with parameters as shown on the graph (solid). The dashed 
lines represent energy shifts caused by the single half-spaces 
of refractive index n\ = 2tt (top) and n s = 2 (bottom) . 



thickness of the slab r is equal to the half-integer multi- 
ple of the wavelength of the dominant atomic transition 
Xji. The existence of resonance effects suggests a physi- 
cal picture of the excited atom as a radiating dipole. The 
resonance and anti-resonance correspond to constructive 
and destructive interference. We have also provided rea- 
sonable approximations in the non-retarded (|75[) and re- 



To provide an additional check on the consistency of 
our calculations we would like to derive equation (|66|) 
by means of electrostatics. We start from the general 
formula derived in [20[ that expresses the electrostatic 
interaction energy of a electric dipole in the presence of 
a dielectric in terms of Green's function of the Laplace 
equation 



AE = l'£{rf)V i V' i G H (r,r') 



(Bl) 



r=r ,r'=r 



Here the sum runs over three components of the dipole 
moment and the subscript H means that only the ho- 
mogeneous correction to the free-space Green's function 
that is caused by the presence of the boundary enters the 
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formula. This ensures that the self-energy of the dipole 
is omitted and guarantees the convergence of the final re- 
sult. The harmonic function G^r(r, r') is a solution of the 
Laplace equation that vanishes for \z\ — > oo. Therefore 
it can be written in the form: 



It has been shown there that, to leading-order, the energy 
shift in the retarded limit can be expressed as 



G h {t,t') 



1 



Ci(k|| 
C 2 (k|| 
C 4 (k|| 



4ire 

,r')e kz 
,v')e kz 
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r')e k * 



+ C 3 (k||,r'> 



z < L/2 
\z\<L/2 ,(B2) 
z > L/2 



ae 1 : 



647r 2 e Z' 



E 



will 



E-'i 



(Cl) 



with k z = <Jk%. + ky. The C coefficients are easily 

worked out by applying the continuity conditions, which 
result from Maxwell's equations, across the interfaces and 
one finds that 



with the coefficients cm j_ given by 



G H (r,r') 
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o Jo 



dkJo(kp)e 



1 



-k(z+z') 



-2kL 



i 2 + 1 n 2 + n 2 

nf-lnf-n s 2 c _ 9fefj 



(B3) 



In 



with p = \/(x — x') 2 + (y — y') 2 . Application of the for- 
mula (|B1[) is straightforward and we easily derive that the 
electrostatic interaction energy of a dipole in a vicinity 
of the layered dielectric is indeed equal to (1551) . 



Appendix C: Retarded limit of the interaction 
energy between an atom and a dielectric half-space 

The interaction between an atom and a non-dispersive 
dielectric half-space has been considered in detail in [9j. 
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part of the integrand about k = 0. For a more rigorous 
treatment see 12111. 



